Variation of fundamental parameters and dark energy. A principal component 

approach 
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We discuss methods based on Principal Component Analysis to constrain the dark energy equation 
of state using a combination of Type la supernovae at low redshift and spectroscopic measurements 
of varying fundamental couplings at higher redshifts. We discuss the performance of this method 
when future better-quality datasets are available, focusing on two forthcoming ESO spectrographs 
- ESPRESSO for the VLT and CODEX for the E-ELT - which include these measurements as a 
key part of their science cases. These can realize the prospect of a detailed characterization of dark 
energy properties almost all the way up to redshift 4. 

PACS numbers: 98.80.-k,98.80. Jk 



I. INTRODUCTION 

Cosmology has recently entered a precision, data- 
driven era. The availability of ever larger, higher-quality 
datasets has led to the so-called concordance model. This 
is a remarkably simple model (with a small number of 
free parameters) which provides a very good fit to the 
existing data. However, there is a price to pay for this 
success: the data suggests that 96% of the contents of the 
universe is in a still unknown form. This is often called 
the dark component of the universe. Whatever this may 
be, all the evidence suggests that it is not composed by 
the protons, neutrons and electrons that we are famil- 
iar with, but it must be in some form never seen in the 
laboratory. 

Current best estimates suggest that this dark compo- 
nent is in fact a combination of two distinct components. 
The first is called dark matter (making about 23% of 
the universe) and it is clustered in large-scale structures 
like galaxies. The second, which has gravitational prop- 
erties very similar to those of the cosmological constant 
first proposed by Einstein, is called dark energy and cur- 
rently dominating the universe, with about 73% of the 
density of the universe 

Understanding what constitutes this dark energy is one 
of the most important problems of modern cosmology. 
In particular, we would like to find out if it is indeed a 
cosmological constant [l|, since there are many possible 
alternatives Q. These alternative models often involve 
scalar fields, an example of which is the Higgs field which 
the LHC is searching for. A further alternative are the 
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so-called modified gravity models (for a review see e.g. 0, 
4]), in which the large-scale behavior of the gravitational 
interaction is different from that predicted by Einstein's 
gravity. 

The main difference between the cosmological constant 
and the models involving scalar fields (which are often 
collectively called dynamical dark energy models) is that 
in the first case the density of dark energy is always con- 
stant (it does not get diluted by the expansion of the 
universe) while in the second one the dark energy density 
does change. One way to distinguish the two possibili- 
ties is to find ways to measure the dark energy density 
at several epochs in the universe. 

Astrophysical measurements of nature's dimensionless 
fundamental coupling constants 0-0] can be used to 
study the properties of dark energy, either by themselves 
or in combination with other cosmological datasets (such 
as Type la supernovas and the cosmic microwave back- 
ground). The concept behind this method is described 
in [9i-tl2| (see also [13| ). It complements other methods 
due to its large redshift lever arm and the fact that these 
measurements can be done from ground-based telescopes. 

Here we revisit this issue and forecast the number 
of well constrained modes of the dark energy equation 
of state parameter using a combination of supernovae 
data and measurements of varying fundamental cou- 
plings at high redshift. We will test and validate the 
reconstruction pipelines by applying them to simulated 
datasets representative of forthcoming high-quality mea- 
surements. This is particularly relevant for the measure- 
ments of varying couplings: the existing spectroscopic 
measurements of the fine-structure constant a I 14| -|l , i 
and the proton to electron mass ratio /z [18l - l20l | typi- 
cally come from observations that were not gathered with 
this purpose in mind, and therefore may be vulnerable 
to considerable uncertainties that are not always easy to 
quantify. 

We will quantify by em ploying Principal Component 
Analysis (PCA) (see e.g. [2l|, [22| ) what improvements 
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will result from the availability of spectrographs like 
ESPRESSO (for the VLT) and CODEX (for the E-ELT), 
which include measurements of a and fi as a key part 
of their science cases. For this purpose we will assume, 
in either case, several scenarios for the datasets of finc- 
structure constant measurements, which will differ in the 
number and precision of the measurements. 

We should stress that we are not proposing new data 
analysis techniques. We are extending PCA methods 
available in the published literature (for type la super- 
novae, lensing and several other contexts in cosmology) 
and studying the feasibility of applying them to a new 
type of datasets (astrophysical measurements of varying 
couplings). 

Briefly, the plan of the paper is as follows. In Sec. 
HTl we briefly review the PCA technique, as it applies to 
our present purposes. In particular we discuss possible 
strategies for choosing the number of components and de- 
scribe how to build the relevant Fisher matrix for the case 
of varying couplings. In Sec. IHII we apply our methods to 
several scenarios relevant for ESPRESSO and CODEX. 
In Sec. IIVI we attempt a reconstruction of the equation 
of state parameter w(z) using a truncation of the high 
frequency modes. Finally in Sec. [V] we present some 
conclusions and highlight future work. 



the value Wi, 

N 

w(z) = J2 w Mz)- (1) 

i=l 

Another way of saying this is that w(z) is expanded in 
the basis 0,, with X = (1,0,0,...), 2 = (0,1,0,...), etc. 

The precision on the measurement of Wi can be inferred 
from the Fisher matrix of the parameters Wi , specifically 
from J (F~ 1 )a, and increases for larger redshift. One 
can however find a basis in which all the parameters are 
uncorrelated. This can be done by simply diagonalizing 
the Fisher matrix such that F — W T AW where A is 
diagonal and the rows of W are the eigenvectors e^(z) or 
the principal components. These define the new basis in 
which the new coefficients an are uncorrelated and now 
we can write 

N 

V)(z) = 22 a i e i( Z ) ■ ( 2 ) 
i=l 

The diagonal elements of A are the eigenvalues Aj (or- 
dered from largest to smallest) and define the variance of 
the new parameters, <7 2 (a,) = 1/Aj. 

A. Building the Fisher matrix 



II. PRINCIPAL COMPONENT ANALYSIS 

PCA is a non-parametric method for constraining the 
dark energy equation of state. In assessing its perfor- 
mance, one should not compare it to parametric meth- 
ods. Indeed, no such comparison is possible (even in 
principle), since the two methods arc addressing differ- 
ent questions. Instead one should compare with another 
non-parametric reconstruction, and for our purposes with 
varying couplings the type la supernovae provide a rele- 
vant comparison. 

A key advantage of PCA techniques is that they al- 
low one to infer which and how many parameters can 
be most accurately determined with a given experiment. 
Instead of assuming a parametrization for the relevant 
observable with a set of parameters born of our theoreti- 
cal prejudices, the PCA method leaves the issue of finding 
the best parametrization to be decided by the data itself. 
This is particular useful in the case of dark energy where, 
apart from the case of a cosmological constant, one would 
be hard pressed to find solid motivations for particular 
parametrizations. 

In Refs. pH . [22l | the PCA approach was applied to the 
use of supernova data to constrain the dark energy equa- 
tion of state, w(z). We start by recalling some of their 
formalism, which we will then generalize for measure- 
ments of the fine-structure constant a. In general one 
can divide the redshift range of the survey into N bins 
such that in bin i the equation of state parameter takes 



We will consider the standard class of models for which 
the variation of the fine-structure constant a is linearly 
proportional to the displacement of a scalar field, and 
further assume that this field is a quintessence type field, 
i.e. responsible for the current acceleration of the Uni- 
verse [23 - i29j . We take the coupling between the scalar 
field and electromagnetism to be 



= ~B F (c(>)F^F^, 



where the gauge kinetic function Bp(<j)) is linear, 



(3) 



(4) 



k = 8irG and C is a constant to be marginalized over. 
This can be seen as the first term of a Taylor expansion, 
and should be a good approximation if the field is slowly 
varying at low redshift. Then, the evolution of a is given 
by 
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For a flat Friedmann-Robertson- Walker Universe with a 
canonical scalar field, cp 2 = (1 + w(z))p ( j > , hence, for a 
given dependence of the equation of state parameter w(z) 
with redshift, the scalar field evolves as 
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where we have chosen the positive root of the solution. 

Let us construct the Likelihood function for a generic 
observable m(zi,Wi,c) — fj,(zi,Wi) + c. For the present 
purposes this can be the apparent magnitude of a super- 
nova, in which case 

/i = 5\og{H d L ) , c = M + 25-51ogif (7) 

or it can be connected to the relative variation of a ob- 
tained with quasar absorption spectra, for which 



H = ln[re(0 - cf>o)} 



c = lnC. 



Then we find 



L(w\M) oc exp 



1 N 

■ E 



mp 



),C % Mm-m F ) 3 



(9) 

where mp means m evaluated at the fiducial values of 
the parameters, mp — mp(z.i,w F ,c F ) and C _1 is the 
inverse of the correlation matrix of the data. 

Defining (3 — c — c F , and integrating the likelihood in 



(8) /3, we obtain the marginalized likelihood 



L(wi) 

where A = and 



exp 



1 N 
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i N 
°y ~ a 2^ 



(10) 
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The Fisher matrix can be obtained by approximating L{wi) as a Gaussian in the theoretical parameters Wi (the 
equation of state in each bin) centered around the fiducial model, and taking the inverse of the resulting correlation 
function. The Fisher matrix turns out to be 
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where the derivatives are evaluated at the fiducial values 
of the parameters. 

We will consider three fiducial forms of the equation of 
state parameter: 

w F {z) = -0.9, (11) 
w F {z) = -0.5 + 0.5tanh(z- 1.5), (12) 
w F {z) = -0.9+ 1.3cxp(-(z- 1.5) 2 /0.1) • (13) 

These three examples cover three classes of the possible 
behavior of w(z). The first corresponds to a constant 
equation of state parameter, the second to a slow tran- 
sition from a dust like component at large redshifts to a 
cosmological constant type at low redshifts and the third 
corresponds to a sharp transition in the value of a scalar 
field occurring around redshift z = 1.5 (see [12J for fur- 
ther discussion). 

We emphasize that the first two parametrizations lead 
to a variation of Aa/a that does not satisfy current geo- 
physical bounds of the Oklo natural reactor [3FJ, HU . (Al- 
though the interpretation of the Oklo bound is not as 
straightforward as one may think, this is not the place to 
discuss it.) Thus for the purposes of the present paper we 
use these 3 parametrizations solely as a proof-of-concept 



representatives of the whole zoo of models, for the pur- 
pose of exploring the reconstruction method. 



III. EXPLORING THE PIPELINE 

We are now in a position to start a forecast analysis of 
the dark energy equation of state for our three fiducial 
forms of w F (z). We take a total number of bins between 
redshift and 4 to be 30. We assume a sample of 3000 
supernovae distributed between redshift and 1.7 (with 
13 bins) with an uncertainty on the magnitude of a m = 
0.11. These numbers are typical of future supernovae 
datasets. For the spectroscopic measurements we use a 
distribution between redshift 0.5 and 4 (with 27 bins) 
and we will consider three different scenarios. Some bins 
overlap so we obtain in total 30 bins. 

We will assume a flat universe, and further simplify the 
analysis by fixing O m = 0.3. This is a standard proce- 
dure, that was followed in the original paper of Huterer 
and Starkman [2ll ] and also in a number of subsequent 
works. The goal of the analysis is to characterize gains in 
sensitivity as future, more precise datasets become avail- 
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able, rather than provide hard numbers. Therefore, al- 
though this is certainly a simplifying assumption, it is a 
legitimate one. This specific choice of Vl m has a negligi- 
ble effect on the main result of the analysis, which is the 
uncertainty in the best determined modes. 

A. Forthcoming datasets 

Based on current plans for ESPRESSO [37| and 
CODEX H (see also Refs. [H H3), we will consider 
three different simulated datasets of fine-structure con- 
stant measurements: 

• A baseline scenario, in which we will assume 30 
systems with cr Aa / a = 6 x 10~ 7 for ESPRESSO and 
100 systems with a Aa/a = 1 x 10~ 7 for CODEX, 
uniformly distributed in the redshift range 0.5 < 
z < 4. This is meant to represent what we can 
confidently expect to achieve in a relatively short 
amount of time once the spectrographs are oper- 
ational (within 3 to 5 years of data acquisition), 
given the current plans for their sensitivity, and it 
will therefore provide the basis for most of our dis- 
cussion. 

• An ideal scenario, in which we will assume 100 
systems with a Aaja = 2 x 10~ 7 for ESPRESSO and 
150 systems with a Aa/a = 3x 10~ 8 for CODEX. 
This is optimistic both in the uncertainty of in- 
dividual measurements and in the number of mea- 
surements. Although several hundred absorbers are 
already known where these measurements can be 
carried out, the sources are quite faint and putting 
together such a dataset would at the very least re- 
quire a very long time. Having said that, our goal 
in considering this case is to obtain an indication 
for the dependence of our results on the uncertainty 
and number of the measurements. 

• A control scenario, which is meant to be some- 
what more realistic from an observational per- 
spective in the sense that we do not assume the 
same uncertainty for all measurements (although 
we still assume that they are uniformly distributed 
in the redshift range 0.5 < z < 4). In this case, 
for ESPRESSO, we assumed the same number of 
sources as in the baseline scenario but with the un- 
certainties drawn from a normal distribution cen- 
tered on <i Aa / a — 6 x 10~ 7 and with standard de- 
viation a Aa / a /2; the same can also be done for 
CODEX. This is a computationally simple way to 
check how the pipeline handles non-uniform uncer- 
tainties and is also a proxy for the effect of redshift 
coverage - an important issue when defining an ob- 
servational strategy. 

We can now use the principal component approach 
to compare the various data sets as probes of dark en- 
ergy. We start by showing the spectra of eigenvalues (or 
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FIG. 1: The error <jj for the five best determined modes for 
the fiducial parametrization (jlll) in the baseline scenario. 
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FIG. 2: The error en for the five best determined modes for 
the fiducial parametrization in the ideal scenario. 



equivalently, the error on the mode i, af = 1/Aj) for 
the different scenarios that we proposed to analyse, for 
parametrization (fTTj) . This is illustrated in Figs. [TJ [2] 
and [3J The shaded region represents a rough threshold 
a = 0.3 below which modes are considered well deter- 
mined, and therefore, informative [34j . By inspecting 
Figs. Q] to [3] it can easily be understood that the base- 
line and control scenarios yield identical results for the 
magnitude of the modes' errors. The ideal scenario, as 
expected, permits that more eigenmodes enter the region 
for which it is considered that modes are informative. In 
order to avoid an overburden of similar figures, we sum- 
marize the eigenvalue spectra for this parametrization 
(fTi) and parametrizations H2]) and (Tj~3]) in Tables HI — [TTT1 
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FIG. 3: The error at for the five best determined modes for 
the fiducial parametrization JTT} in the control scenario. 
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TABLE I: The table indicates, for the fiducial parametrization 
JTT} , how many modes have an error below the threshold value 
a = 0.3 considering the baseline, ideal and control scenarios. 



It can be noticed that the baseline and control scenarios 
for ESPRESSO cannot bring further information when 
compared with supernovae alone. Only in the ideal sce- 
nario and in combination with supernovae, can the num- 
ber of useful modes raise above 2. The only exception 
is for parametrization (fTTj) where the ideal ESPRESSO 
has 4 good modes. CODEX, however, will offer a set 
of high precision data which can on its own provide a 
large number of informative eigenmodes in either of the 
three scenarios. We can see that the results are fairly in- 
dependent of the fiducial parametrization chosen though 
parametrization (jllj) allows, in general, a higher number 
of modes considered informative for a given combination 
of data. 

Let us now look at the evolution of the eigenmodes for 
parametrization and compare with the modes using 
supernovae data alone. Because only the first two modes 
are considered informative, a reconstruction of the equa- 
tion of state parameter using only supernovae can only 
be made reliably up to redshift 0.5 using the principal 
component approach. This value corresponds roughly to 
the position of the second mode's peak in Fig. SJ We 
are interested in evaluating whether this limit can be re- 
laxed with the help of QSO data using the ESPRESSO 
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TABLE II: The table indicates, for the fiducial parametriza- 
tion (|12|) , how many modes have an error below the threshold 
value a — 0.3 considering the baseline, ideal and control sce- 
narios. 
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TABLE III: The table indicates, for the fiducial parametriza- 
tion (|13[) . how many modes have an error below the threshold 
value a — 0.3 considering the baseline, ideal and control sce- 
narios. 



and CODEX spectrographs. For the baseline scenario 
we see that only CODEX will permit a large number of 
modes considered useful. Inspecting Fig. [5] we see that a 
reconstruction of the equation of state parameter could 
be in principle possible up to redshift 2 corresponding 
roughly to the third peak of the third mode. The conclu- 
sions are similar for the ideal and control scenarios, see 
Figs. [5] and [7J In particular, only for the ideal scenario 
would ESPRESSO add information on the nature of dark 
energy. The form of the fiducial function w F (z) is not 
determinant in the results and we would have reached 
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FIG. 4: The four best determined eigenmodes for 
parametrization using only supernovae. Solid line, 

dashed line, dash-dotted line and dotted line correspond to 
first, second, third and fourth modes, respectively. 
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FIG. 5: The four best determined eigenmodes using 
parametrization (j 1 1 p for the baseline scenario. Solid line, 
dashed line, dash-dotted line and dotted line correspond to 
first, second, third and fourth modes, respectively. 



FIG. 7: The four best determined eigenmodes using 
parametrization fTTJ for the control scenario. Solid line, 
dashed line, dash-dotted line and dotted line correspond to 
first, second, third and fourth modes, respectively. 
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FIG. 6: The four best determined eigenmodes using 
parametrization (llip for the ideal scenario. Solid line, dashed 
line, dash-dotted line and dotted line correspond to first, sec- 
ond, third and fourth modes, respectively. 



IV. RECONSTRUCTING THE EQUATION OF 
STATE PARAMETER w(z) 

Folowing Refs. [2l], [22j one can now attempt a recon- 
struction w{z) by keeping only the most accurately de- 
termined modes (the ones with largest eigenvalues). To 
do this, we need to decide how many components to keep. 
We must point out that the weak point of this procedure 
consists in neglecting the high frequency modes. In a 
more recent analysis [34j and perhaps more robust ap- 
proach, all modes are kept and a correlation function 
describing fluctuations from a fiducial model is chosen. 
This method allows a more accurate reconstruction of 
the dark energy parameter over the whole range of red- 
shifts covered by observations. We decided, however, to 
follow Refs. [U, [22| given the simplicity of the methods 
described there. This will be sufficient to support our 
main conclusion i.e., that a combination of supernovae 
and quasar data will improve on testing the time depen- 
dence of dark energy. 



A. Selection of components: risk vs. normalization 



the same conclusions had we illustrated the eigenmode 
evolution of parametrizations (fT2|) or (fT3| . We must em- 
phasize, however, that the choice of threshold does play a 
role. Had we chosen a slightly higher value, say a = 0.5 
then the baseline scenario of ESPRESSO could indeed 
bring new information on the nature of dark energy when 
comparing to supernovae observations. 



One may argue that the optimal value of modes M to 
be kept corresponds to the value that minimizes the risk, 
defined as [2l| 

risk = bias 2 + variance, (14) 

with 

N 

Uas 2 {M) = ~ w F {zi)) 2 , (15) 

i=l 
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FIG. 8: Reconstruction of the equation of state parameter 
JTT} using only supernovae with the minimization of the risk 
method. 

where the notation w means that the sum in ([2J runs 
from 1 to M, and 

N M 

variance — J]] J]] °" 2 ( a j) e j( z i)- (16) 
i=i j=i 

The bias measures how much the reconstructed equation 
of state, w TCC (z), differs from the true one by neglecting 
the high and noisy modes and therefore, typically de- 
creases as we increase M. The variance of w(z), however, 
will increase as we increase M, since we will be including 
modes that are less accurately determined. 

An alternative way to decide on the number of optimal 
modes is to choose the largest value for which the error 
is below unity, or equivalently, the RMS fluctuations of 
the equation of state parameter in such mode are 

((l + w(z)f)=a?<l. (17) 

Having thus determined the optimal number of modes, 
we proceed with the normalization of the error follow- 
ing Ref. [35| such that a 2 = 1 for the worse determined 
mode and normalize the error on the remaining modes 
by taking 

a\ ai ) ^ aKa,) = f ffi (18) 

1 + <7-{Oli) 

The PCA allows us to optimize an experiment towards 
the range in redshift we are interested in. In our work, we 
will use a combination of supernovae and quasar absorp- 
tion lines to understand how well the equation of state pa- 
rameter w(z) will be constrained with forthcoming data 
on cosmological variation of fundamental parameters ob- 
tained with the spectrographs ESPRESSO for the VLT 
and CODEX for the E-ELT. 

In Fig. [5] we illustrate the reconstruction using only 
supernovae for our fiducial model (|TTj) and in Figs. [5] and 



[10] we show the result of the reconstruction for the two 
spectrographs and the two methods for the selection of 
the number of components. We can observe that, for this 
fixed number of bins, the reconstruction obtained using 
supernovae is only accurate up to redshift s» 1. In par- 
ticular, because we neglect the poorly determined modes 
which are the ones with high amplitudes for bins of large 
redshift, the reconstructed equation of state parameter 
tends to zero for large redshift. This unavoidable feature 
of the PCA truncation method can be confused with a 
real increase in the e quat ion of state at high redshift. As 
it was pointed out in [2l|, had we reconstructed 1 + w(z) 
instead, then w{z) would approach —1 for large redshifts. 
Ideally we would like to extend the survey to large red- 
shifts but unfortunately we only expect data from super- 
novas up to redshift z w 1.7. 

Measurements from the quasar absorption lines, which 
are available for a larger redshift interval, provide in 
general a more reliable reconstruction. For our fiducial 
parametrizations of w(z), these datasets can give a quali- 
tatively accurate account of the evolution of the equation 
of state parameter to fairly high redshifts. 

Comparing the various fiducial models for the same ob- 
servational dataset shows that (as one would expect) the 
redshift up to which the reconstruction remains accurate 
depends in part on the correct underlying model, specif- 
ically on whether its equation of state remains close to a 
cosmological constant or approaches a dust-like behavior. 
However, comparing the CODEX and ESPRESSO cases 
show that one can go deeper in redshift by increasing the 
sensitivity of the measurements, since that allows one to 
add components to the reconstruction. In particular, the 
truncation problem mentioned earlier becomes less prob- 
lematic as the reconstructed equation of state parameter 
no longer approaches zero for large redshifts. This is a 
true statement also for the alternative parametrizations. 

The combination of supernovae with quasar absorp- 
tion lines data further improves the determination of the 
equation of state parameter. In particular, we can now 
obtain information on w(z) all the way from z up 
to z rj 3. The reconstruction using CODEX, benefiting 
from an almost one order of magnitude improvement in 
the sensitivity of the QSO data points, is substantially 
better than the one obtained with ESPRESSO. 

We can also compare the two methods of determining 
the optimal value of modes to keep in the reconstruc- 
tion. We have seen before that the minimization of the 
risk method is a compromise between having an accu- 
rate equation of state reconstruction and having a small 
error bar in this reconstruction. The normalization of 
the error on the modes method, however, makes use of 
our prior prejudice that variations of w(z) larger than 
unity are unlikely. We observe, from comparing Fig. [5] 
to Fig. [Till to which correspond different methods, that 
the latter method picks more modes, which leads to a 
more accurate reconstruction. Since we are including ad- 
ditional modes with progressively larger errors, the recon- 
structed equation of state in this case also has larger error 
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FIG. 9: Reconstruction of the equation of state parameter 
JTTJ in the baseline scenario with the minimization of the risk 
method. 
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FIG. 10: Reconstruction of the equation of state parameter 
(|ll|l in the baseline scenario with the normalization of the 
error on the modes method. 



bars. In other words, the normalization method provides 
a more conservative and accurate approach, while the 
risk method provides (appropriately) a more aggressive 
approach. 

From a practical point of view, the advantage of PCA 
techniques is that they may provide us with a simple com- 
putational tool to continually optimize an observational 
plan. In other words, given an ongoing observational 
campaign in which one has already observed a certain 
number of sources (with given uncertainties) one can use 
these tools to simulate improved datasets with the goal 
of determining how to best spend the remaining tele- 
scope time (reducing the uncertainty in measurements of 
particular sources) so as to achieve the best possible con- 
straints on these models, given any relevant observational 



limitations. 



V. CONCLUSIONS 

In the first part of our study we considered the thresh- 
old value presented in Ref. [34] as a mean to determine 
how many modes can be considered good and informative 
in describing dark energy. This value, a < 0.3, is some- 
how arbitrary but serves our purpose to indicate that fu- 
ture precise measurements of the fine-structure constant, 
specially with CODEX, can add information on the na- 
ture of dark energy up to high redshifts when comparing 
to supernovae measurements. 

In the second part of this work we have compared re- 
constructions of the equation of state parameter w(z) of 
dark energy using a principal component analysis. To 
this effect we used a combination of expected supernovae 
data and quasar absorption measurements of the fine- 
structure constant a, expected to be available with the 
spectrographs ESPRESSO and CODEX. We considered 
several possible datasets and also two methods of choos- 
ing the best determined modes of the principal decompo- 
sition and studied the effect of the size of the error bars 
on the reconstruction. 

Our analysis indicates that the normalization of the 
error on the modes method appears to give more accu- 
rate (closer to the fiducial value) but less precise (more 
conservative errors) reconstructions with respect to the 
risk minimization procedure. We also conclude that a re- 
construction using quasar absorption lines is expected to 
be more accurate than using supernovae data. However, 
since the two types of measurements probe different (but 
overlapping) redshift ranges, combining them leads to a 
more complete picture of the evolution of the equation of 
state parameter between redshift zero and four. 

A natural extension of this work is to include in addi- 
tion cosmological measurements of ix = m p /m e , the ratio 
of the proton and electron masses [20| . Measurements of 
fi will be fewer than those of a but will all be at high 
redshift and fairly precise and are, therefore, expected to 
reduce the errors on the reconstruction of w(z). 

Although in this work we have focused on combining 
two datasets, we should also point out that they can also 
be used separately to provide independent reconstruc- 
tions. Comparing the two reconstructions will then pro- 
vide a consistency test, specifically for the assumption 
on the coupling between the scalar field and electromag- 
netism (given by Eq. 0J. If so one can also obtain a 
measurement for the coupling parameter £. A more de- 
tailed treatment of this well as the application of 
the method to existing datasets, is left for forthcoming 
work. 

Since a PCA reconstruction involves a truncation, it 
will not be precise at high redshifts. How much of a 
problem this is depends on the behaviour of the true 
equation of state: the reconstructed equation will always 
approach at high redshifts, while the true one may or 
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may not do so in the redshift range being probed by 
the detaset under consideration. A relevant question for 
PCA studies is therefore how deep in redshift can one 
confidently go. In most previous works one has a fiducial 
model for w{z) that approaches by around z ~ 1, in 
which case supernovas perform well; however this is an 
optimistic assumption that need not be true. One of our 
key points is that, not knowing a priori what the correct 
equation of state will look like, there is strong interest in 
trying to reconstruct as deep as possible in redshift, and 
varying couplings are a possible way of doing that. 

One must bear in mind that one should not compare 
the parametric-free PCA constraints with those obtained 
using parametrized equations of state. Perhaps a more 
direct comparison is with redshift binned w(z). A good 
benchmark can be found in Fig. 17 of the work by Kowal- 
sky et al. [36J, where one sees how poorly the high red- 
shift behaviour of the equation of state is measured by 
current data. 

Finally, our main conclusion is that, when one com- 
pares like with like, the inclusion of varying constants 
data allows a reliable reconstruction to be carried out 
to significantly higher redshifts. Admittedly, for these 
methods to be competitive will require very good qual- 
ity measurements, which do not presently exist (although 
they are expected to be available in a few years), but this 



point applies both for the supernovas and for the varying 
couplings. In any case, the method is interesting on its 
own right, and should be further studied as the science 
cases of future facilities is developed. 
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